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1. INTRODUCTION 


In this paper we consider two types of phenomena which a- 
rise in the theory of low energy fission due to the finite ve- 
locity of a fission process: diabatic mode coupling (DMC) for 
collective coordinates and energy loss of the fission mode via 
diasipative processes: If the fission motion were infinitely 
slow one could describe the wave function of a fissioning nuc- 
leus within the zero order Born-Oppenheimer approximation (BOA ), 
i.e. the wave function would separate into the simple product 
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у (A,, ) = U (a) W( q,»;«) a 


Here х is the fission coordinate, q stands for the collective 
degrees of freedom of the fissioning nucleus and x for the re- 
maining internal coordinates. U ( ) is the probability amplitude 
for finding the fission coordinate between & and а da, while 
И (х, а; < ) is the wave function in the space of all other coordi- 
nates; it depends parametrically on х , In the case of BOA the 
wave function F (x, q: c) belongs to the lowest eigenvalue of the 
internal Hamiltonian being compatible with the prescribed con- 
served quantum numbers, The first step towards the inclusion of 
effects of the finite fission velocity consists in expanding the 
wave function in а basis which allows for excited states in the 
X,q-space and for subsequent modifications of the fission wave 
function itself: 


y (a, . x) = ý C nk U mk (A) W- (q, <) (2) 


m, K 
Here Na (x,q$ ) is an eigenfunction of the internal Hamiltonian 
Hint (x, q; &) which depends parametrically on : 


Н, 609.59) Wm (Galia) = би () Wm (q, x ja) (3) 


The fission wave function depends on the quantum numbers m of 
the internal state via the energy Eu and on the quantum num= 
bers К which characterize the fission mode (energy, angular mo- 
mentum etc.). The coupling between the excited modes 4 шк (& . qX ) 
results from the operator 14%) which is obtained from the part 
of the kinetic energy operator which acts on the parametric de- 
pendence of Ta on & , Since this parametric dependence of LA 
on ais directly related to the deviation from the simple adia- 
batic wave function of Eq.(1), we call the coupling leading to 
the form of Eq.(2) diabatic mode coupling. 

The IMC would adequately describe the fission process, if 
in Eq.(2) we could really perform the sums over all relevant 
quantum numbers m, i.e. sums over quantum numbers related to 


Nonadiabatic and Dissipative Effects ... 363 


collective and intrinsic coordinates. Since this is mathematical- 
ly unfeasable ore can divide the problem into two distinct rarts: 
In the first step one ignores the intrinsic degrees of freedom 
and performs a DMC calculation in q-space™, In the second step 
one takes account of the coupling between the intrinsic coordina- 
tes т and the collective coordinates < and д. The most  promi- 
nent effect of this coupling is the thermalization of the energy 
of the collective motion. Two mechanisms can be distinguished: 
(i) The coupling between the intrinsic coordinates x and the 
collective coordinates q leads via intrinsic excitations to a 
damping of the collective mode (e.g. the asymmetry vibration). 
This is the same mechanism that leads to the damping of excited 
collective states of fission-stable nuclei, as for example in 
£iant resonances, This damping depends on the fission motion in- 
directly via the available energy and the energy dependence of 
the damping width. (ii) The coupling between the intrinsic coordi- 
nates x and the fission coordinate X leads to a direct excita- 
tion of intrinsic degrees of freedom which in turn results in a 
direct dissipation of the fission energy. The dependence on the 
fission mode comes via the velocity of the fission process and via 
the fission energy which determines the phase space for intrinsic 
excitationse 

In our treatment we ignore the change of nuclear structure 
parameters due to the heating up of the fissioning nucleus, though 
this will be probably important for the inertia parameters, These 
problems will be considered in further investigations. 

Furthermore our treatment is limited to small fission ener- 
gies, otherwise diabatic level crossing would become an important 
effect which upsets the basis of our description. 


2. BASIC DEFINITIONS AND NOTATIONS 


In recent years many calculations of potential surfaces and 
inertia tensors for fission configurations of nuclei have been 


Xin the numerical application we deal with only one collective co- 
ordinate belonging to the mass asymmetry oscillation. In principle 
more collective coordinates could be included. 
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performed (1-5) and methods have been developed which allow the 

determination of adiabatic paths to fission [5]. To describe the 
fission process in the adiabatic approximation one generally in- 
troduces a number of collective coordinates: The fission coordi- 
nate & related to the mean elongation, a coordinate des descri- 
bing a possible mass asymmetry, a necking coordinate Oder and 

eventually further coordinates depending on the degree of sophi- 


stication one is aiming at. 
In this paper we use tne collective coordinates of ref. [5]; 


1.2. 


a (a, + C, + С,) /2Ко, 
9 * (a,b - a,b; la, bit arb) s (AA: ДА,+А,). 


With these coordinates the classical energy of a system of par- 
ticles becomes 


44 5 vi 93i E 708), (2.1) 


where 171 and (a) denote the mass tensor and the potential 
energy. After quantization the Hamiltonian takes the form 


¥ -h * 1 ` 
H=-7 2 |M| 2 nj M e Z vs), (2.2) 
Ы Sunt | 


where tne determinant of the mass tensor is |M| = det м; 4) and 
м, ut = х, 
19 c3 е 
In the two dimensional space with the coordinates &, q 
the mass tensor is 
ka ( Maa Mag 
M. = 
tj \ Mga Maq / 


(2.3) 


We chose а coordinate system 9; = (g) in which the nass ten- 
зог is diagonal На = Ka = ©. 
Ye denote the diagonal elements of the mass tensor 


Nonadiabatic and Dissipative Effects s.s 365 


(2.4) 
gd Mg. 
They depend on « and qe 
То describe the fission process we have to solve the eigen- 
problem with the collective Haniltonian 


3 r 2. 2 Ме г. | (2.5) 
2 Mag d I Ma 39 | Mq 24 : 
The variables & and q describe the fission and asymmetry modes 
respectively. The masses М. («,q) and Ko (N,) аз well as the 
potential (,d) are obtained from the asymmetric two center 
shell model [5]. 

The scalar product of the eigenfunctions of the Hamiltonian 
H is defined with the metric 


D(a.q) -N (, q) Ha, J ; (4,/ч.) «f [9 (9,9) D(x,q)du dq 


To obtain a metric which does not depend on masses and coordinates 
and gives orthogonal functions 


(9,/9,) -& 9, (a. q) de dg 


we perform now an unitary transformation of the Hamiltonian and 
its eigenfunctions у (,) with the help of the function | (a :3) 


9 (, q) = yD(a3) 9 (x, q), 
A (a, q) VDA N, q) үре. 


Then the Hamiltonian takes the form 


(2.6) 


accent 2 L 2. | (vq + Vg (6,9). (2.7) 


dą Ma да 29 Mą д9. 


Neglecting the small scalar term Yo (= „d) (Appendix А) we split 
up the Hamiltonian in the following way: 
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= fe- На , (2.8) 
where 
a uy rdi "o 
nh. alii oc E € | (2.8a) 
» 2 
n (2.8) 


Ha = 2 dą Ма(е,9) 99 


In the adiabatic approximation one assumes that M depends only 
weakly on the fission cordinate &. Ye therefore define an adia- 
batic Hamiltonian for the internal system 


ad т. : 
Ho (arg) 4 a TORT * V(q;«). (2.9) 


where М (q) is an average mass for the  g-motion. наё defines 


the eigenfunctions of the q-mode 


4% (458) Wn (G32) = бт (=) Мм (q; =) 


In order to define the adiabatic part of 7, we introduce the 
average mass 


Ма (a) / g ;) Ма (2,9) W(q;«)/|Wlq;a)|* dą 


Then gat is defined by 


„ 
< 2 da Mal) aa’ (2.10) 


The eigenfunctions of the total adiabatic Hamiltonian 


(2.11) 


н“ (aq) стад 4 на“ 


can be written аз: 


CET = Це (x) w, (gj (2.12) 
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The adiabatic Schrödinger equation is 


ad Р f : 
(rat нае) ujęta) м, (9:9) EU es) eons 
which yields the following equation for the fission wave function: 
ad 
(TĘ + Ce) Une (a) * EU; (a) (2.14) 
The energies f(a)» £ (а), 8 3 KC define the fission po- 


tential when the fission takes place with the internal sy st en EM 
in the states n 1,253, 


i 
i 
1 
i 
t 
= 
a бб Św 
Fig, 1. Effective potential Ё, (a) for 
different internal states 81% n 


3. DESCRIPTION OF THE MODEL . 


Without dissipation the Hamiltonian in Eq. (2) is equal to 
H = H + qRAd, The energies En (a) are а sum of two terms: 
E 1209) + J (a) where J () represents the fission potential 
and Tp ( ) m = 0,1,2, is the energy of the m-th asymmetry pho= 
non. Then we include dissipation the asymmetry mode is no longer 
stable because it can now thermalize, In principle this thermali- 
zation could de treated by the explicit introduction of all the 
degrees of freedom into which the phonon can decay. However, in 
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practice this procedure is quite unfZeasable. The sell known re- 
cipe to get around this problem is FesLbach's effective operator 
method (10]: One continues to work in the space defined by the 
unperturbed phonons (Р-зрасе) and taxes account of the remaining 
space into which the phonons decay (2-space) by the introduction 
of an effective energy dependent Hamiltonian. In this way the 
phonon energies become complex quantities and lead ultimately to 
a complex fission potential which is, of course, the analog of 
the complex optical potential for nucleon-nucleus scattering: in 
the course of the fission motion the amplitudes | m,E > must 
be dazred because they are partly scattered into more complicated 
States. 

70 consider this picture in more detail we write for the 


operator H in Eq. (1.3) 


Я. + (g. a) = V (a) e (g, x, c) 


and replace the phonon operator by an effective operator PH, E (E Р 
which is defined in P-space 


РН, (E)P = РН, (x) P+ РҮ (E-QHQ) QvP 3.1) 


Here Hola) is the unperturbed Hamiltonian defined by 


Но () Wa +5, (%) Wa 


and V is the interaction which couples the phonon to more com- 

plicated configurations», There are two distinct contributions to 

this coupling: 

(i) For fixed « the finite velocity of the asymmetry oscillation 
leads to a coupling V. = V(q,x) between the coordinate q 
and the internal coordinates x; if the asymmetry motion is 
approximated by an octupole vibration the form of this cou- 
pling is well known from the physics of low lying collective 
states (see e.g. ref. [11]. It leads to the usual damping of 
the collective motion as it is seen in another context in 
the width of giant resonances, 


(11) Since the velocity of the fission motion is finite, there 
will be also a direct coupling Ух = V(0,x) between the 
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fission coordinate а and the internal degrees of freedom x. 
The damping associated with this type of coupling has some 
similiarity to friction because it is directly related to 
the velocity of the fission motion. 


The determination of the precise form of the coupling terms 
Va and M^ and the evaluation of Herr?) is а aifficult pro- 
blem which is far from being solved, For our purposed where the 
effective interaction enters only in an average way the following 
procedure appears to be reasonable: 

The operators Ух and Ya are sums Of single particle ore- 
rators. We apply them to the wavefunction W, (X, , 4) and obtain 


ub 
ui W„(q,x;a) = v Wry, pe g. e à A „ (5.2) 


where Lm у (49% %) is the wavefunction of а state where either 
a ph-pair u, y has been created "on top of" the phonon state Я 
or of a noncollective ph state «which is orthogonal to the col- 

lective state Wa Taking matrix elements of H f S) we obtain 


m 


ŚW, | Нед | Wy gy n X LE He | Wa Da? 


| i 

P <w; | Mal Wi iy T ŚW, uy quate | Pn us М, 
ży 

g> <W. "Ми: y Mum | eee: "T Уи; uv | \ | Wy № 
n (5.5) 


The matrix elements in this equation are defined by integrations 
over x and q, they depend parametrically on c. 

In the following we shall consider only diagonal elements 
of the effective operator. Clearly in principle there vill be 
also nondiagonal matrix elementa, so that one would have to per- 
form a diagonalization of Berg C) in the space | Ww. > + Due to 
phase canceilation effects the nondiagonal matrix elements should 
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be on the average much smaller than the diagonal ones. Therefore 
the main effect of the second an third term on the r.h.s. of Eq. 
(3.3) will be on the diagonal matrix elements. Also for phase con- 
cellation reasons we have neglected the interference terms between 


Va and У. 
In order to evaluate the matrix elements of the resolvent 


(E-QHQ we proceed as follows: 
G) Fe mt CW [o цу "бди бум. 


which is familiar from shell model calculations. 


(ii) We replace the sum over д and y by an integration over 
the corresponding energy variable introducing the density of 
states 9 (w). Че can then write Eq. (5.5) in the form 


2 s 2 
| Маана | + IVa mia) | 
Cni (ED n a (nie [dmg (je) Е-Е; (a) -w e i8 | 
6.4) 
Va (1, , © ) 3 and var а) 2 are averages of the squares of 


the matrix elements. 
For the real and imaginary part of (3.4) we obtain 


ginaf (ра) Va (i zna] 


О te — 77777 


(3.58) 


Jm Św; | Hay (ENW: X =- T g (Е -£&);«) [I (н) Мн] 
(5. 5b) 


The real part implies a renormalization of the phonon energies 
which can be discarded for the following. The imaginary part 
makes the phonon energies complex or, in other words, leads to 
the appearence of an imaginary part in the fission potential; 
1% describes the weakening of the amplitude 9:8 (<) 11 (dx, c) 


in the course of the fission process. 
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It is a difficult problem to estimate the coupling matrix 
elements Y, and Та. Te propose the following approxirate 
treatment: For fixed coordinate а and the constituent par- 
ticles see a mean potential Ua,» 905 fn) acting on particle n. 
A displacement aAa Ar a, wą leads to a change of the 


au au J 
potential (a= ao) ją а= A and (q - a) 24 q = 9 


respectively. 
Рог the asymmetry coordinate q the perturbation 


(9 = 95) ui | 9 = 9% leads to the conventional form of the 


particle-phonon coupling and the resulting damping width will Бе 
the width of the asymmetry phonon. 

For small deformations, i.e. for small values of а, the 
quantity | CE can be estimated rather reliably from the widths 
of giant resonances. The matrix elements are typically of the 
order of 0.1 Мет. For large deformations there ray be considerable 
deviations from these values; at the moment ме do not consider 
such deviations but they should be taken into account in future 
more refined studies. 

The change of the fission coordinate & leads to dissipation 
of the fission kinetic energy via onebody and twobody dissiration. 
The resulting friction force F = ~ 1% is related to the icasi- 
nary part of the optical potential via the danping time: The 
friction force leads to a damping of the energy of the fission 
motion with the relaxation time t, = m/2T where m is the cor- 
responding inertia parameter. On the other hand an imaginary part 
* in the optical potential implies a damping of the energy with 
the relaxation time tz = A/Z. The counter part of the friction 
force in the optical potential is therefore an imaginary part 


w /m à (3.6) 


For low energy fission the motion is slow enough that the 
friction coefficient can be calculated in the framework of linear 
response theory as discussed in refs.[12 $ 15]. А calculation of 
T in this spirit has been performed in ref.[14] for moderate de- 
formations (around the second minimum of the fission potential) 
and extended to large deformations by Marcev (priv. communication). 
The calculated friction coefficient is an energy dependent quan- 
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tity r (м). The large width of the energy distribution - of the 
order of 10 He- implies a memory time of the order of 10 "22; 
which is short compared to characteristic times o? the fission 
motion. Therefore memory effects can be safely neglected and 

T (w— 9) Gefines the friction coefficient in the ordinary sense, 


4. RESULTS AND DISCUSSIONS 


4.1 In order to separate the effects of mode-coupling and of dis- 
sipation we have first made a calculation where dissipation 
was completely neglected. The calculations are done for the 
nucleus 2569. The ground state fission barrier is double 
humped, in the region of the first maximum and the second 
minimum the system follows a symmetric path and nonadiabatic 
effects are believed to be small before the second maximum. 
For this reason we assume the interaction operator quad to 
vanish to the left of the turning point атут (see Fig. 1) 
belonging to ground state fission. The aave function Y 2 we 
are looking for contains in addition to the adiabatic wave 


function ф 92 c*her contributions boś»: 


ad s LĄ т! „a - * T 
4, (6.9) (257 fate SMS UT) m 
121 


The coefficients с, (Z,E'; are determined from a variational 
principle in the following way: пе make the ansatz 


бие) 35 a gh (Ее FERES] 42 


k 2 4 U. 


which allows the system to gc off shell in the fission energy 
up to |E - E'| = 4.5 the function h, is the Rermite poly- 
nomial of order К. Since we go up to Ках = 10 the form(4.2) 
garantees a great flexibility for the functional dependence 
of e, on the off-shell energy E’ ~ Е. The coefficients ak 
are determined from the requirement 
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a * 2 2 
day LU (н-Е) Уст (a) da dą | »Q (4.3) 
which leads to the following system of linear equations for 
the coefficients any? 
Kin N r Р 
$ oy A aud uus Ue (4-4) 
К =4 121 , 


The normalisation factor n' (&) = | vag (< |^? is introduced 
for the following reason: In the classically forbidden region 
the'fission wave function decays very rapidly vith increasing 
A. Therefore, without the factor | U(a )| 2, the integral in 
eq. (4.3) receives practically no contribution from regions 
of X which are far away from the turning point х1 11 and the 
variational method is consequently only sensitive to details 
of E wavefunction in the vicinity of qa 111 · Multiplication 
by n? (A) corrects for this defect and leads to equal weights 
of the internal wavefunction in the entire intervall of the 
X integration. 

The matrix A is determined as 


MT : [da 4% J (a.) дих (9) * (9) 
sith 
Fen 4) * (H- e) cf iI K, :£ ep S; 


ad 
where H-E acts on q gg (& . q) · 


The vector b is given by 
cfe dq 9* (q) q (2,9) vf (a) 


The width a5 is determined from the requirement that the 
norm of H B) Wg(«,q) should be as small as possible. 
We found 45 ~ 0,6 MeV. 
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Ammon RE ˙ —˙—V AO = 


The probability distributions for the masses of the fission 
fragments are obtained from the solution of eq. (4.1) for 
a & a ac! 


[4 
CACIK = |o (a q (4.5) 


Unfortunately the scission point & = 2.5 can not be reached 
in our calculation. But we found that beyond & ~ 1.75 the 
mass distributions do not change significantly (see fig. 2); 
we therefore calculated the mass distribution from eq. (4.5) 
with & = 1.85. 


Rig. 2. Dependence of 
mass distribution on 
the fission coordinate 


for 25 a 


In order to give an idea of the level spacings of the phonon 
states associated with the asymmetry mode we show in fig. 4. 
the energies of the 10 lowest phonons as a function ofa . 

It is seen that they follow without much fluctuations the po- 
tential surface belonging to the groundstate. It is clear 
from the figure that the number of phonons that one has to 
include will increase almost linearly with the fission energy. 


In fig. 3 we show for the nucleus 2360 the mass distribution 
probability following from the adiabatic wave function - which 
ig of course independent of the fission energy - and the mass 
distribution probability following from nonadiabatic mode 
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Fig. 3. Dependence 
of mass distribution 
on the fission co- 


ordinate r 256; 


coupling for E = 2, 4 and 6 MeV. 
It is seen that with increasing 
fission energy the curves become 
broader and develope a shoulder 
towards more symmetric fissior. 
The small bump for A ~ 133 is 
prcbably a numerical effect 
being due to a still too small 
phonon space. 

The conclusion from these 
findings is that nonadiabatic 
coupling is unimportant for 
Е <3 MeV but becomes quite 
important for higher fission 


Fig. 4. Effective fission poten- energies. However; for ener- 
tial for the 10 — symmetry y 
phonon states. gies above 4 MeV our approxi- 


mation begins to fail because 
diabatic level crossing becames а very important effect which 
dominates the process of energy dissipation [15]. 


4.2 A second type of calculations was done in order to isolate 
the effects of the dissipative terms in eq. (2.8). They were 
done only in the region a > Ay since the increasing fission 
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velocity in this region is believed to lead to large friction 
effects. In addition, the classically forbidden region is quite 
uninteresting in this context since the damping of the fission 
wave function due to barrier penetration is at any rate much 
stronger than the damping due to the imaginary part of the fis- 
sion potential. In these calculations we have neglected the 
nonadiabatic coupling pnad, therefore there is only one channel 
- in which the q-mode is in its ground state - which is treated 
explicitely. А 

This procedure is justified а posterori by the findings 
that: a) DUC is only slightly affected by the frictional датріг 
of the x -mode. b) The loss by Die is small compared to the 
friction induced energy loss of the fission motion. 

In ref. [14] the friction coefficient was calculated with 
a coupling paraneter Ea = 0.03 меу"! which was fitted in ref. 
[16] to optical potential data in an energy range of 5-0 Hey. 
However the main contributions to y comes from much smaller 
energies for which the value Ге = 0.03 leads to notoriously 
too small spreading width of the quasiparticle states, Че there- 
fore have chosen Г= 0,15 which yields much better results for 
the spreading width of states close to the Fermi surface. Ji ch 
cz = 0.15 мет”! яа AE Tir % 50 A. Рог the mass para- 
meter we took m = 100 Hey which corresponds to 70 a. n. u. 

It ig interesting to compare these figures to the results 
of pure onebody dissipation obtained from the mall formula [17]. 
Taking a mean value for the multipolarity I = 5 - which ap- 
pears reasonable for the large deformation between saddle and 
scission point - one obtain Ty, & 300% and m z 80 &^/MeY. 

The value of Tir = 50 h comprises the effects of the 
moving walls and of residual interactions, і.е, one- and two- 
body dissipation. Therefore it might look surprising that one- 
body dissipation within the wall formula yields a much larger 
friction coefficient. However, one must realize that the wall 
formula is based on a single particle picture for the motions 
of the nucleons where all nucleons participate in the same 
manner in the interaction with the moving walls, In reality 
the single particle like states are quasiparticle states. They 
possess only a finite lifetine oT oop. „Only for states 
close to the Fermi energy is this lifetime large compared to 
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a typical nucieon passing tines of 10-22 sec. For states lying 
5-10 MeV below the Fermi energy the imaginary part of the auasi- 
particle self-energy attains ~ 5 MeV [16] corresponding to 

ora. p. = 10722 sec. Рог increasing distance from tne Ferzi 
energy ata. P. falls rapidly belom the passing time, For the 
interaction process underlying the wail formula this means that 
during the time intervall in which a quasiparticle enters into 
the diffuse surface region of the nucleus and leaves it again it 
dissolves into a multitude of long range (n + 1) particle -n hole 
configurations - which partly correspond to dynamical degress of 
freedom of the nuclear surface - zhile the single particle con= 
ponent dies out very quickly. Therefore the effective nuuber of 
particles which have to be counted in the application of the wall 
must be considerably smaller than the nucleon number A. A quanti- 
tative evaluation of the effect of the quasiparticle decay on one 
.onebody dissipation is difficult but qualitatively the discre- 
pancy between Tir and Tag is not astonishing. 

With our choice of parameters we obtain YT = 0.8 MeV for 
the imaginary part of the optical potential. 

Due to the large collective mass the fission wavefunction 
can be calculated in the "KB approximation. The changes which 
become necessary due to the presence of an imaginary part in the 
potential are rather straightforward. Че discuss them in App. C. 

The value of Я = 0.8 Ме7 implies a strong damping of the 
elastic channel. For instance, for ground state fission the pro- 
bability damping factor (see egs. (C5) and (C6)) 


c 
' 
2 " 1 a da е 
p exp [ 57 |ą( )| (4.6) 
dy 
becomes Poamp z 6.1074; consequently there are practically no 
cold fiscion fragments left when the system arrives at the scis- 
sion configuration. 
The energy absorption is governed by the relation 


«|? E (a) exp [- Ia], 
pit 


E (a+8a) = E(x) exp [- 2—3 
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where va is the collective velocity. Vith Е = 5 v? оре 


obtains 


= z œ T Va (4.8 ) 


which is the classical relation following from the Rayleigh dis- 
sipation function. Our optical model treatment together with a 
quasiclassical treatment of the motion leads us back to the 
classical problem of the motion of a particle which moves in an 
external field and feels viscous damping. The quantum nature of 
the problem enters only via the determination of the friction 
coefficient, 

The numerical treatment of the fission motion leads to the 
following results (for ground state fission ): 


22 
(i) The time required to go from Жуу to * 18 g 22.10 ae 
вес, 
(ii) The gain in kinetic energy of the fission mode in going 
from “ту to agi is „ 10 Key. It is seen that the motion 
is not as crrepy as predicted by the wall formula. 


(iii) The intervall fron ту to X + 1,85 where DMC is most ef- 
fective is passed within 7.10722 sec. The collective velo- 
city at © = 1.85 is K = vy = 0,039 «1022 sec compared 
to vq 0. 050. 1032 вес”! for | = 0. The collective ener- 
БУ logs at «= 1.85 is ~ 5.5 MeV. The mall relative change 
which results from switching on friction shows that the 
damping of the fission mode has only a small effect on DMC, 


(iv) The typical period of the asymmetry oscillations Ta 
15 1077? sec is only slightly smaller than Tyy » so that 
there is just enough time for me period while the system 
goes from Ary to Ase’ During this tíme the asymmetry 
phonons do not have sufficient time to fully thermalize 
because with a typical width of 0,1 to 0,2 MeV - estimated 
from speading widths of low lying collective states - the 
thermalization time lies between 30 and 60 10722 вес. 
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APPENDIX А 


Estimation of the geodetic potential 


On the Hamiltonian defined in Eq. (2.5) we perform а 
transformation with the help of the function 


YD =УМа (%9) Mq (9,9). (к.1) 
іп the following way: 
-4 
Hy 8 — (0 H(V0) VO g:EY5 4 (4.2) 


- -4 
The new Hamiltonian and wavefunctions are Н = VD H (vo) 


and ф = (D Ф respectively. Apparently the transformation con- 
serves the energy eigenvalues. The transformation of the first 
term 


Бан ы . 
Ta * Ма № Ka Ma д“ в 


a 
r 
Mm 


gives 


2 { 9 i 
FR (he 4) 
where 
a aca Em 4 з. 4 Z tvo] 
Va- M, UD at Vo VD (ôa Fo da AM 
z 5 4 JS NS 02. 
From the second term "4°" мам ay ET 29 (A. 6) 
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FF 
one obtains 7 2 ^ Wa gehe 4.7) 


LM 


VS . © ża 5715 o 


żę у 


Finally one gets for the transformed Hamiltonian 


vieil GM: ea als nem ae vl oe 
H= 2 da Ma dua 39 Mq dq "Aven oe — 


The quantity V,(a,q)= UM. е v£) is the last term appear- 
ing in Eq. (2.6) on the r.h.s. Numerical estimates show that it 
can be neglected compared to the other terme in H, 
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APPENDIX B 


„valuasion of the fission wave function in WKB approximation 
fhe collective Hamiltonian connected with the coordinate 


jescribing the fission mode is assumed in the following form: 


? 2 
| Y 1 не) $e 5 LUla) - jugi v Pata 


where М(х) is the collective mass parameter and U(a ) is the 
collective potential which consists of the two parts: 


Uta) =М (м) + i W(a) (B.2) 
The real part V of the potential is the fission barrier and the 
imaginary part Я 1s responsible for nuciear dissipation. 


The first order VKB approximation gives the following general 
solution for U (a): 


U(a)* С quart et fs a Jda] exp (4 f pida- Zorc tg E 
(B.5) 
er Nc fg / er eo fp nnt n] 


as will be shown in Appendix C. 


Here 
M(a) [E-V(a) + c- Mah z И? («) ] (B.5a) 
and 
E M («Y W(a) 
3 p (e? . 30) 


The collective velocity v is given by 
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—— ————— DDD 


| р(х) | 
У (а) = "R * (a) 


Let us make the following assumption about the real part of 
collective potential 


(8.5с) 


const=V (d zr) for a < Agr 
(G) =- < Va(*) for Я gr $ & C4 (8.4) 


constaV(a,.) for «5» its 


where a sr corresponds to the ground state configuration and 
0 ас to the scission point, The potential (4) is shown in Pig. 1. 
The last assumption in (4) (fora > & ue. means that the omit 
the Coulomb tail. 

The imaginary part of the potential is assumed in the fol- 
lowing form 


0 for © (m, 
* () < (а) for A LA 4 Wee (B.5) 
о for X > Ase 


It means that we assume the pfesence of nuclear dissipation in 

the region from the right turning point a, (see Fig. 1) to the 
scission configuration A c° n (a) and W„(a ) are the micro- 
scopic estimates of the real and imaginary parts of the collec- 
tive potential respectively. This choice of the potential ensures 
that the incoming (fora < X gr) and outgoing (а > Qa вс) wave 
functions are plane waves. The WKB wave function which corresponda 
$0 the potential (4) and (5) ia given by (see e.g. Landau and 
Lifshitz, Quantum Mechanics) 
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a 

$ ofi вв.) tor a <a, 
ac 

| а 
u 1 1 | 6 
(ау = ye p| El f pda for a <A C, (B.6 ) 
| ac 
* 


Ka- 
Ww 


. jest | 


for ауа; 
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APPENDIX С 


WKB for complex potential and variable mass 


Ге start from the Schrödinger equation 


[5 = em du Че! 3) u (a) = U (6.1) 
with the complex potential 
U (a) = V(a)+i W(a) 
Taking the aneatz 
u (a) ~ exp (16 (a)/ fi) (С.2а) 
* 6 (&) 6 (a) +Еб"( а) . 25) 


one obtains the following equation for 6 (a) 


dé Y? ih Gu d?6 _ ие: | F 


а) лба] 2Mia) OEG 


2M 


(C.3) 
Making the usual expansion 


в Y eo 
B" 6^, 2 6. ) RÓ 


we obtain after a lengthy but straightforward calculation the 
following expressions for 6 ©) ana 6 


6 (a) =t f pla) аа t C, 


6a) =: [q(a)da + Co (0.4) 
(9' ых Ve) + аа) 

6 (a) ? [п ie) © + C5 

6" (a) =- 5 осм Mr IS >. "С. 
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where Cy Co» cz; e. are arbitrary integration constants 
and р(х) and q(«) are defined as 


P (a) e [E-vi«) ДЕ-Уе nea) | 


M (a) W (a) 
9 (a) = - pla) 


дуја one obtains the following form for 6 (д) = 6 ©) ta ) + 
665 
a) 


(С.5.) 


Е: * are tq 365 ec tC, 


Ll 
+i E ahd. ms + Co-hC4 


M (=) 


This gives the WKB wave function 


Ula) +) at | È Ner- fatis) | pa- Зас ta nj 
+ D,exp Aq exp[- € Греда 2: arctq " 


where D, and D, are arbitrary constants which heve to be chosen 
in accordance with the boundary conditions. In our case where we 
have а plane incident wave coming from the left we put D, z0 
(recall that W< 0). 
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(v) Comparing the effects coming from DMC and damped fission 


1. 


motion one sees that the energy loss is primarily determined 
by the frictional type of damping of the fission mode; com- 
paratively little energy goes into the excitation of asym- 
metry phonons. On the other hand the latter process governs 
the fission fragment distribution. It is formed in an early 
stage of the decent towards scission where the motion is 
slow and the energy loss via friction is small so that the 
frictional energy dissipation has little influence on mass 
distributions. 

Everyghing that has been said is valid only for small 
fission energies probably not higher than 3 or 4 MeV, Above 
this limit nonadiabatic level crossing becomes so important 
that it invalidates the whole approach so that a completely 
different theory must be envisaged. 
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STRESZCZENIE 


Praca dotyczy oceny rozkładu mas jąder fragmentów rozszcze- 
pienia i dyssypacji energii jądra 2560 Obliczenia ва wykonane 
w przybliżeniu adiabatycznym, ale dyskutowane są także: a) po- 
prawki do rozkładu шав ze względu na nieadiabatyczne sprzężenie 
z fononami asymetrycznych wibracji, b) tłumienie ruchu do roz- 
szczepienia poprzez dyssypację jedno- i dwuciałową ze wspórczyn- 
nikami tarcia obliczonymi na drodze mikroskopowej. 
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РЕЗЮМЕ 


Работа посвящена оценке распределения масс осколков де- 
ления и диссирации энергии деляющегося ядра U 36, Расчёт ос- 
нован Ha адиабатическим приближении но мы учитываем тоже: 

а) добавки к распределении масс из за неадиабатического спа- 
ривания фононов ассиметрических вибраций; б) задерхание ne- 

ления из за 1- и 2-частичной диссипации, где параметр трения 
определён по микроскопической теории. 


